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Abstract 

We study the nonclassicality of photon-added squeezed vacuum (PASV) and its decoherence 
in thermal environment in terms of the sub-Poissonian statistics and the negativity of Wigner 
function (WF). By converting the PASV to a squeezed Hermite polynomial excitation state, we 
derive a compact expression for the normalization factor of m-PASV, which is an m-order Legendre 
polynomial of squeezing parameter r. We also derive the explicit expression of WF of m-PASV and 
find the negative region of WF in phase space. We show that there is an upper bound value of r for 
this state to exhibit sub-Poissonian statistics increasing as m increases. Then we derive the explicit 
analytical expression of time evolution of WF of m-PASV in the thermal channel and discuss the 
loss of nonclassicality using the negativity of WF. The threshold value of decay time is presented 
for the single PASV. 

PACS number(s): 03.67.-a, 03.65. Ud, 42. 50.Dv 



1 Introduction 

Nonclassical states play an important role in quantum computation and quantum imformation pro- 
cession 111 . The nonclassicality of quantum states can be well-described by some nonclassical prop- 
erties, such as sub-Poissonian photon ststistics [2], squeezing in one of the quadratures of the field 
0, and negativity of Wigner function (WF) [4]. Especially, the partial negativity of WF implies the 
highly nonclassical properties of quantum states and is often used to describe the decoherence of 
quantum states. Many experimental schemes have been proposed to generate nonclassical states 
of optical field. Among them, subtracting photons form and/or adding photons to quantum states 
have been paid much attention because these fields exhibit an abundant of nonclassical properties 
and may give access to a complete engineering of quantum states and to fundamental quantum 
phenomena J5j [6j |7J [8j [TOl 1111 [121 113B . For example, a single photon-subtraction squeezed vacuum 
(PSSV) has been experimentally prepared with a pulsed and continuous wave squeezed vacuum 
Bill 1121 IT3H . It is very similar to quantum superpositions of coherent states with small amplitudes 
H141 PT5H . As another example, a single-photon addition was experimentally performed by Zavatta et 
al [6], which unveils the nonclassical features associated with the excitation of a classical coherent 
field by a single light quantum. For the single PSSV, its nonclassical properties and decoherence was 
investigated theoretically in two different decoherent channels (amplitude decay and phase damp- 
ing) by Biswas and Agarwal H1611 . They indicated that the WF losses its non-Gaussian nature and 
becomes Gaussian at long times in amplitude decay case. 
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On the other hand, the combination of the photon addition and subtraction has been successfully 
demonstrated in Ref. fl8j |9j ■ In Ref.|[8j|, photon addition and subtraction experimentally have been 
employed to probe quantum commutation rules by Parigi et al. In fact, they have implemented sim- 
ple alternated sequences of photon creation (addition) and annihilation (subtraction) on a thermal 
field and observed the noncommutativity of the creation and annihilation operators. It is interest- 
ing to notice that subtracting or adding one photon from/to pure squeezed vacuums can generate 
the same output state, i.e., squeezed single-photon state H17L However, the resulting states ob- 
tained by successive photon subtractions or additions are different from each other. For instance, 
successive two-photon additions [a^ 2 ] and successive two-photon subtractions [a 2 ] will result in the 
same state produced by using subtraction-addition (a^a) and addition-subtraction (aa f ), (also see 
section 2 below) respectively. As far as we know, the nonclassicality and decoherence of arbitrary 
number photon-added squeezed vacuum states (PASV) in a dissipative channel has not been derived 
analytically in the literature before. 

In this paper, we shall investigate the nonclassical properties and decoherence of single-mode 
PASV which is optically produced single-mode non-Gaussian states. This work is arranged as follows. 
In section 2, we introduce the single-mode PASV and discuss its nonclassicality in terms of sub- 
Poissonian statistics and the negativity of its Wigner function (WF). By converting the PASV to a 
squeezed Hermite polynomial excitation state, we derive a compact expression for the normalization 
factor of PASV, which is an m-order Legendre polynomial of the squeezing parameter r, where m 
is the number of added photons; and then derive the explicit analytical expression of WF for any 
photon-added number m and find the negative region of WF in phase space. We also show that 
there is an upper bound value of r for this state to exhibit sub-Poissonian statistics which increases 
as m increases. In section 3, we derive the explicit analytical expression of time evolution of WF 
of the arbitrary PASV in the thermal channel and discuss the loss of nonclassicality in reference of 
the negativity of WF. The threshold value of the decay time corresponding to the transition of the 
WF from partial negative to completely positive definite is obtained at the center of the phase space, 
which is independent of the squeezing parameter. We show that the WF for single PASV has always 
negative value for all parameters r if the decay time nt < \ In 2 ^+ 2 (see Eq.Q40D below), where 
n denotes the average thermal photon number in the environment with dissipative coefficient n. 
Conclusions are involved in the last section. 

2 Single-mode Photon added squeezed vacuum state 

Various photon states have been generated by the micromaser and WFs of some cavity fields can be 
measured by a scheme based on interaction between cavity fields and atoms H181 |l9l . As described 
in Ref. [12011 when an excited atom passes through a cavity field which is in a squeezed vacuum state 
then their interaction may produce photon addition (excitation) on the squeezed vacuum state — the 
excited squeezed vacuum state. 

Therectically the single-mode PASV can be obtained by repeatedly operating the photon creation 
operator on a squeezed vacuum state S (r) |0), i.e., 

\r,m)=N r>m ^ m S(r)\0), (1) 

where |0) is single mode vacuum, and N r . m is the normalization constant to be determined, a) is the 
Bose creation operator, and S (r) is the single-mode squeezing operator S (A) = exp[i (ra* 2 — ra 2 )] 
H211 122H with r being the squeezing parameter. 

2.1 Single-mode PASV as the squeezed Hermite polynomial excitation state 

Under the transform of S (A) we see (r) S (r) = a) coshr + asinhr, thus we can reform Eq.© 
as 

\r,m)=N r , m S{r) (r) a^ m S(r) |0) 

= N r ^ m S (r) (cr cosh?' + a sinhr) m |0). (2) 



On the other hand, using the operator identity [23 
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where H m (x) is m-order single variable Hermite polynomial whose definition is 
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which indicates that the PASV is equivalent to a squeezed Hermite-excited vacuum state. 
Further using the generating function of H m (x), 



H m (x) = — exp (2xt - t 2 ) 
and Eq.©, the normalization factor N r _ rn can be derived by 
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where in the last step we have used the Baker-Hausdorff formula e^ a e va ' < = e" a1 e^e^. Using the 
newly found generating function of Legendre polynomial H24l I25H (see Appendix A), 
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we can derive the compact expression for N r , m , i.e., 

2V-£ = ml cosh m rP ro (cosh r) , (8) 

which is just the result in Ref.[20] derived by the mathematical induction method. In particular, 
when m = 1, (Hi (x) = 2x) i.e., the single-photon added case, we see that |r, 1) = S (r) |0) is just 
a squeezed single-photon state; on the other hand, for the single-photon subtracted case H171I26H , 
the state is aS (r) |0) = S (r) (r) aS (r) |0) = 5 (r) (acoshr + a+ sinhr) |0) -> 5 (r) a f |0) , which 
indicates that adding a single-photon to the squeezed state has the same impact as annihilating 
a photon from the squeezed state. While for m ^ 2, the case is not true (also see Fig.l). For 
example, successive two-photon additions [a^ 2 ] and successive two-photon subtractions [a 2 ] will 
result in the same state produced by using subtraction-addition (a^a) and addition-subtraction (aa)), 
respectively, i.e., a) 2 S (r) |0) -► a)aS (r) |0) , a 2 S (r) |0) -► aa f 5 (r) |0) . In Ref.ll27ll, two PSSV is used 
to generate the squeezed superposition of coherent states with high fidelities and large amplitudes. 



Combining Eqs.flT]) and ([8]) we can conveniently calculate the average photon number a^a in 
PASV, 

(r, m\ a) a \r, m) — (r, m\ aa) |r, to) — 1 
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where we denote £ = cosh r for simplicity, and 

(r, m| a^ 2 a 2 |r, to) = (r, to| (a 2 a^ 2 — 4aa^ + 2) |r, to) 
TV 2 iV 2 



A/2 ju2 
iv r,m+2 Jv r,m+1 



thus the Mandel's Q-parameter can be obtained by substituting Eqs.® and dlQD into Q = 
(a^a). In particular, for single-photon-added case m = 1, Eqs.© and fllQD reduce to 
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(r,l|a t2 a 2 |r,l) = — ^, (12) 



thus the Q-parameter with m — 1 is given by 

Q = 3slnh22r -1. (13) 
3 cosh 2r — 1 

From Eq. (ll3D we find that Q becomes negative for to = 1 which is satisfied for the squeezing 
parameter r < 0.46 similar to the result of Q for single-photon subtracted case 111611 . In order to 
see clearly the variation of Q-parameter with r, we show the plots of Q-parameter in Fig.l, from 
which one can clearly see that Q-parameter becomes negetive (m ^ 0) when r is less than a certain 
threshold value which increases as m increases; while for to = 0, Q is always positive. This implies 
that the nonclassicality is enhanced by adding photon to squeezed state. We should emphasize that 
the WF has negative region for all r, and thus the PASV is nonclassical. In our following work, we 
pay attention to the (ideal) PASV in a thermal channel. 

2.2 Wigner function of PASV 

In order to discuss the decoherence properties of PASV in thermal environment, in this subsection 
we shall derive the analytical expression of Wigner function for PASVS. For single-mode case, the 
Wigner operator is defined as ||28ll 
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where \q) is the coordinate representation, Q\q) = q\q). Thus the Wigner function of PASVS |r, to) 
can be calculated by W (q,p) = (r, m\ A (q,p) \r, to) . Using Eq.Q, we can see 



W = sinh '" 2r (°l H m (-*\l^ a ) St to A S W H m L^^aA |0) . (15) 



Figure 1: (Color online) The Q-parameter as the function of squeezing parameter r for different m = 
0,1,2,3,4,29,30. 



On the other hand, noticing that the single-mode squeezing operator (r) has its natural expres- 
sion in coordinate representation H29H , i.e., 5* (r) = e r l 2 f^° oD dq\qe r ) (q\ , leading to (r) \q) = 
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= A(qe- r ,pe r ) . (16) 

Then substituting Eq. (ll6D into Eq. (I15D and using Eq.© as well as the coherent state representation 
of Wigner operator, 
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where a = (q + ip)/^/2 and \z) = exp(za^ — z*a) |0) is the coherent state H301 1311 . we can put 
Eq. (115D into the following form, 
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where 



(3 = — i/3y2 cothr, /3 = (ge _r + ipe r ) /y/2 = a cosh r — a* sinh r, 
and in the last step, we have used the integration formula 

/2. 
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In order to further simplify Eq.( 118D , expanding the exponential item e 2tTCOthr as sum series and 
using Eq.© we have 
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Noticing the recurrence relations of H m (x), 
then the Wigner function of |r, to) is given by 
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where /? and (3 are shown in Eq.fll9D . Eq. (l23D seems a new result (not reported in the literature be- 
fore), related to single-variable Hermite polynomials. In particular, when the photon-added number 
to = 0, 1 and noticing that P (coshr) = 1 and Pi (coshr) = coshr, Eq.( |23D reduce to, respectively, 
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Eq.( 124D is just the WF of squeezed vacuum state, a Gaussian in phase space; while Eq. (l25D corre- 
sponds to a non-Gaussian WF in phase space due to the presence of non-Gaussian item 2(q 2 e~ 2r + 
p 2 e 2r ) — 1. It is clear that the function W m =i (q, p) becomes negative in phase space, when 2(q 2 e~ 2r + 
p 2 e 2r ) < 1. 

Using Eq. (l23D we show the plots of WF in the phase space in Figs. 2 for different squeezing 
parameters r and photon-added numbers to. One can see clearly that there is some negative region 
of the WF in the phase space which implies the nonclassicality of this state. In addition, the squeezing 
effect in one of the quadratures is clear in the plots (see Figs. 2a and Figs. 2b), which is another 
evidence of the nonclassicality of this state. The WF has its minimum value for m = 1,3 at the 
center of phase space (q = p = 0) (see Fig.2(a) and (d)). The case is not true for m = 2 (see Fig2.c). 



3 Decoherence of PASV in thermal environment 
3.1 Model of Decoherence 

In this section, we consider how this single-mode state ([T]) evolves at the presence of thermal envi- 
ronment. In thermal channel, the evolution of the density matrix for the m-PASV can be described 
by El 

= k (n + 1) (2apa) — a)ap — pa) a) + nn (2a) pa — aa) p — paa)) , (26) 

where k represents the dissipative coefficient and n denotes the average thermal photon number of 
the environment. When n = 0, Eq.( |26t reduces to the master equation describing the photon-loss 
channel. The corresponding time evolution density operator can be obtained as 

oo 

Pr, m (*) = e Kt+r ° 2 M km M\ v (27) 



Figure 2: (Color online) The Wigner functions of PASV for (a) m = l,r = 0.3,(b) m = l,r = 0.8, (c) 
m = 2, r — 0.3 and (d) m = 3, r = 0.3. 



where po = \r, m) (r, m\ is the initial density matrix; M k ,i and Ml ; are Hermite conjugated operators 
(Kraus operator) with each other, 
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and T + , T_ and T are determined by 
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It is not difficult to prove the M k j obeys the normalization condition J^Ti^o ^1 i-^k,l = 1 by using 
the technique of integration within an ordered products of operators 



3.2 Wigner function of the PASV in a thermal channel 

The evolution formula of WF of the PASV can be derived as follows [13511 
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where W (a, a*,0) is the WF of the initial state. Eq. d30D is just the evolution formula of WF in 
thermal channel. Thus the WF at any time can be obtained by performing the integration when the 
initial WF is known. 

In a similar way to deriving Eq.( 123D , substituting Eq. d23D into Eq.( 130D and using the generating 
function of single-variable Hermite polynomials (O and Eq. (l20D , we finally obtain (see appendix B) 
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where we have set 
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E = — {AD - 2D* sinh 2r) , (35) 
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Eq. (13 ID is just the analytical expression of WF for the PASV in thermal channel. It is obvious that the 
WF loss its Gaussian property due to the presence of single-variable Hermite polynomials. In partic- 
ular, at the initial time (t = 0), noting \/C (2n + 1) T ->• 2, G -> 1, F/G -> 0, and E/VG -> ft = 
-iy/2 coth r(C cosh r - (* sinhr) , as well as j$ \B\ 2 - ( ^f 1)T -> -2 |C| 2 cosh 2r, (S* 2 + B 2 ) -> 

(C 2 + C* 2 ) sinh2r, Eq.(f3lT) just does reduce to Eq.([23]), i.e., the WF of the PASV. On the other hand, 
when Kt -> oo, noticing that T -> 1, S -> 0, C -> 4, D -> 0, S/VG 0, G -> 1, F -> 4cothr,as 

well as H m (0) = -jrS m _23, then Eq. (l3ll becomes ff ( 2 ^ +1 ) e"^ 1 , which is independent of 

photon-addition number m and corresponds to the WF of thermal state with mean thermal photon 
number n. This indicates that the system state reduces to a thermal state after an enough long time 
interaction with the environment. 

In addition, for the case of m = 0, single-mode squeezed vacuum, Eq. Q31D just becomes (Hq(x) — 

1) 

W m=0 (C,e,t)=me-^ 2+ ^* 2+ ^, (38) 

where m = ,Vc(L + i } t is the normalization factor, D = j^TW - Serb ^ gffl ' E( l-^ 
denotes a Gaussian distribution function — the WF of single-mode squeezed vacuum in the thermal 
channel; while for m = 1, single-photon added case, its WF in the thermal channel is given by 
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In Fig.3, the WFs of the PASV with m = 1 are depicted in phase space with r = 0.3 and n = 1 
for several different k£. It is easy to see that the negative region of WF gradually diminishes as the 
time Kt increases. Actually, from Eq. (|33D one can see that C > 0, so at the center of the phase space 
{a~ a* =0), when F < 2 cothr leading to the following condition: 

1, 2n + 2 , AM 
nt < Kt c = - In . (40) 

2 2n + l' 

which is independent of the squeezing parameter r, there always exist negative region for WF in 
phase space and the WF of PASV is always positive in the whole phase space when nt exceeds the 
threshold value nt c . 

In Figs. 4 and 5, we plot the variation of WF in phase space for different n and r, respectively. It 
is found that the partial negativity of WF decreases gradually as n (or r) increases for a given time. 
The squeezing effect in one of the quadrature is shown in Fig. 5. In addition, for the case of large 
squeezing value r, the single-photon added squeezed state becomes similar to a Schodinger cat state 
(see Fig. 6). The WF becomes Gaussian with the time evolution. In principle, by using the explicit 
expression of WF in Eq. (l31D . we can draw the WF distribution for any photon-added case in phase 
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Figure 3: (Color online) The Wigner functions of single-photon-added squeezed vacuum states in phase space 
for r = 0.3, n = 1 at (a) nt = 0.05, (b) nt = 0.1, (c) «i = 0.2 and (d) Kt = 0.5. 

space. For instance, for m = 2, there are two negative regions of the WF, which differs from the case 
of single PASV (see Fig. 7). The absolute value of the negative minimum of the WF decreases as nt 
increases, which leads to the complete absence of partial negative region. 

4 Conclusions 

The nonclassical properties and decoherence of single-mode PASV in a thermal environment have 
been investigated. A compact expression for the normalization factor of PASV is derived by convert- 
ing the PASV to a squeezed Hermite polynomial excitation state. It is shown that the normalization 
factor is just an m-order Legendre polynomial of the squeezing parameter r. We also derived the 
explicit analytical expression of WF for any photon-added number m and found the negative region 
of WF in phase space. We also show that there is an upper bound value of r for this state to exhibit 
sub-Poissonian statistics which increases as m increases. Then we considered the effects of deco- 
herence to the nonclassicality of PASV when interacting with thermal environment. For arbitrary 
number PASV, we derived the explicit analytical expression of time evolution of WF and presented 
the loss of nonclassicality in reference of the negativity of WF. The threshold value of the decay 
time corresponding to the transition of the WF from partial negative to completely positive definite 
is obtained. We find that the WF has always negative value for all parameters r if the decay time 
< i In §§±f for single PASV. 
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Figure 4: (Color online) The Wigner functions of single-photon-added squeezed vacuum states in phase space 
for r = 0.3 and Kt = 0.05 with (a) n = 0, (6) n = 1, (c) u = 2, (d) rl = 5. 




Figure 5: (Color online) The Wigner functions of single-photon-subtracted squeezed vacuum states in phase 
space for n = 0.1 and nt = 0.05 with (a) r = 0.03, (6) r = 0.5, (c) r = 0.8, (d) r = 1.5.?? 
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Figure 6: (Color online) Wigner function of single photon-added squeezed vacuum states in phase space for 
r = 0.8, n = 0: (a) Kt = 0.1, (b) Kt = 0.2, (c) = 0.3, and (d) nt = 0.7. 



Figure 7: (Color online) The Wigner functions of photon-added squeezed vacuum states in phase space for 
r = 0.7, n = 1, m = 2: (a) Kt = 0.1, (6) nt = 0.2. 



Appendix A: Derivation of Eq.© 

Recalling that the newly found expression of Legendre polynomial 
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though it is different in form from Eq.(Al), they actually are equal. For instance, what we list in the 
following equations (the left) from Eq.(Al) are equal to what we list from (A2) (the right) 
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H 










x 2 ) 


H 


3 
2 




x 2 ) 



and 



m = 4, P 4 (x) = x 4 
m = 5, P5 (a;) = a; 5 



1 + 5 1 



or 
1 



1 + 3 1- — +_ 1- _ 



15 



1 - 



3 2 1 

-a; ; 

2 2' 

5 3 3 

2 X ~ 2 X ' 



= - (35a; 4 - 30a; 2 + 3) 



= - (63a 5 - 70a; 3 + 15a 
8 



(A3) 



(A4) 



We emphasize that the new form in Eq.(Al) cannot be directly obtained by some series summation 
rearrangement technique from the orginal definition (A2) . 

On the other hand, noting the following relation 

Q2m 

^—exp(-t 2 -T 2 + 2xTt)\ tT=0 

\n+l 



E 

i,Lfe=0 



(— k d 2m 



n 



(2x) 



[m/2] 



2 m m \ V 

^-^ 02n 



dt m dr m 
m \ 



„m—2n 



t,T=0 



" Q 2 2n (nly (m - 2n)\ 



(A5) 



and then comparing Eq.(A5) with Eq.(Al) we can obtain Eq.©. 
APPENDIX B: Derivation of Wigner function Eq.([3l]) of PASV 

In this appendix, we present the details for deriving the Wigner function Eq.( 131D . Substitut- 
ing Eq. (l23D into Eq. (l30D and noticing Eq.J19D as well as the generating function of single-variable 
Hermite polynomials (O, we have 



w(c,e,t) = Nj2 



(-2 cothr)' f d 2 a _ 2 \± 



^I![(m-Z)!] 



pfyr 2 1 ol cosh t — ct * sinh r | 2 



\H m -i(—i (a cosh r — a* sinhr) v2 coth? 



1 f)2m-2l 



am' ™ (-2 cothr)' & 

^ll[(m-l)\] 2 dv m - l dT m ~ l 



f d 2 a f ,2 

/ exp | -A \a\ + Bia + B 2 a* 



(a 2 + a* 2 



) sinh 2r j 



r—v—O 



where A and B are given by Eqs. fl34D - (l35D . 

2 sinh" 1 r 



N = 



7T (2n + 1) T 2 m P m (cosh r) ' 



and 



B\ = B* — i2 (v cosh r + t sinh r) V2 cothr, 
Bi = B + i2(v sinh r + r cosh r) V2 cothr. 

Further using the following integral identity H37H : 

J ^ exp (C \z\ 2 + £z + r)z* + fz 2 + , 9 z* 2 ) 



(Bl) 



(B2) 



(B3) 
(B4) 



VC 2 - 4/5 



exp 



C 2 - 4/ 5 



(B5) 



whose convergent condition is Re(C ± / ± g) < 0, and Ref ^±j±| J < 0, we can put Eq.(Bl) into the 
following form: 



w(c,C,t) 



N 



-2 cothr) d 



l f)2m-2l 



^ (2n+1)T ho 11 t( TO " O'] 2 dv m - l dT m ~ l 



- «p \-v 2 - t 2 + ^,B X B 2 + -^{B\ + Bl) sinh2r 1 



(B6) 



where C and I? are given by J34l)-(l35l). and 

B X B 2 = B* B + 2 (vD + tD* ) + 8 (r 2 + v 2 ) cosh 2 r + 8tv cosh 2r coth r, (B7) 
B\+B\= B* 2 + B 2 - 32ru cosh 2 r - 8 (u 2 + r 2 ) cosh 2r coth r - AtD - AvD* , 

then substituting Eqs.(B7)-(B8) into Eq.(B6) yields 



(B8) 



, T r/> ,* x N , 2 KI 2 A| H |2 , gmt2r f B «a , R 2\ y^-y ( — 2 COth r) 

W (C,C , *) = ^= e ' 2 "+ 1 ) T e^ |iJ| +— n— ^ +a ■> \ - 



^![(m-/)r 



d 



2m -21 



g v m-ig T m 



- exp {-G (v 2 + t 2 ) + 2vE + 2tE* + tvF} t=v=q . 



(B9) 



1 o 



where E, F, and G are given by Eqs. J35D - Q37D . In order to further simplify Eq.(B9), expanding the 
exponential item e TvF as sum series and using Eqs.© and the formula 



dv 
we see 



V e -Gv 2 +2vE 



= G m ' 2 H m (e/Vg) , (BIO) 



w (C) C \ t) = JL e -Mw+v\B\ 2 +^{B^+B>) j2 (- 2cothr ^ 

Vc l=0 i\[{m-iy.} 2 

- rpk E)2fc p>2m-2l , „ 

x £ £ £ e -G(« 2 +r 2 )+2^£;+2rE* 

^ fc! d(2E) k d{2E* f dv m - l dT m -> 
= N c _ 7 J^ + A lBl 2 +slI ^ (B ,2 +B 2 ) ™ (-2cothr)'G m -' 

^ «[(m-0!] 2 



T=V=0 



^ pfc n2fc , N ,2 

xV- r j- H m -i[E/VG)\ . (Bll) 

^ *! d(2E) k d(2E*) k m l \ l )\ 

Then using Eq.(f22T) yields Eq.([3TT). 
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